Numerical Solution of the Mode-Coupling Equations for the Kardar-Parisi-Zhang 

Equation in One Dimension 
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We have studied the Kardar-Parisi-Zhang equation in the strong coupling regime in the mode- 
coupling approximation. We solved numerically in dimension d = 1 for the correlation function 
at wavevector k. At large times t we found the predicted stretched exponential decay consistent 
with our previous saddle point analysis in [Phys. Rev. E 63, 057103 (2001)], but we also observed 
that the decay to zero occurred in an unexpected oscillatory way. We have compared the results 
from mode-coupling for the scaling functions with the recent exact results from Prahofer and Spohn 
|cond-mat/0101200| for d = 1 who also find an oscillatory decay to zero. 



PACS numbers: 05.40.-a, 64.60. Ht, 05.70Ln, 68.35.Fx 

The Kardar-Parisi-Zhang (KPZ) Q equation is one 
of the most important phenomenological equations in 
physics. It is essentially a non-linear Langevin equation 
and was proposed in 1986 as a coarse grained description 
of a growing interface. It is the simplest generalization 
of the diffusion equation which includes a relevant non- 
linear term, and probably as a consequence of this the 
KPZ equation also arises in connection with many other 
important physical problems (the Burgers equation for 
one-dimensional turbulence J2j, directed polymers in a 
random medium [^[jH) etc.). 

The KPZ equation in the context of a growing inter- 
face describes it by a single valued height function h(x, t) 
on a d-dimensional substrate x 6 5R d is: 



d t h(x, t) = vV 2 h + ^(Vh) 2 + r/(x, t) . (1) 

The first term on the right of Eq. (1) represents the 
forces which tend to smooth the interface, the second 
describes the non-linear growth locally normal to the 
surface, and the last is a noise term which mimics the 
stochastic nature of the growth process [pi, usually cho- 
sen to be Gaussian, with zero mean and second mo- 
ment (?7(x, t)r)(x.', t')) = 2DS d (x. - x')5(t - t'). The equal 
time interface profile is usually described in terms of the 
roughness: w — y/ (h 2 (x, t)) — {h(x, t)) 2 which for a sys- 
tem of size L behaves like L x f(t/L z ), where f(x) — » 
const as x — ■+ oo and f(x) ~ x x ^ z as x — * 0, so that 
w grows with time like t x ' z until it saturates to L x when 
t ~ L z . x and z are the roughness and dynamic exponent 
respectively. 

Above two dimensions, there are two distinct types 
of solution of the KPZ equation. In the weak coupling 
regime (A < A c ) the non-linear term is irrelevant and 
the behavior is governed by the Gaussian (A = 0) fixed 
point and z = 2. The strong coupling regime (A > A c ), 
where the non-linearity is relevant (and A c = for all 
d < 2) is characterised by exponents which are not known 
exactly in general dimension d. ^,From the Galilean in- 



variance (invariance of Eq. ([!]) under an infinites- 
imal tilting of the surface) one can derive the relation 
X + z = 2, which leaves just one independent expo- 
nent. For the special case when d = 1, the existence of 
a fluctuation-dissipation theorem gives the exact results 
X = 1/2, z = 3/2. 

Almost by definition there is no small parameter for a 
systematic perturbative treatment of the strong coupling 
regime. One is forced either into numerical studies of the 
KPZ equation, which are naturally difficult for dimen- 
sion d greater than two, or into "ad hoc" approximations. 
The best known of these is the so-called mode-coupling 
approximation 0-^], in which in the diagrammatic ex- 
pansion for the correlation and response function only 
diagrams which do not renormalize the three point ver- 
tex A are retained. One of the purposes of this paper is 
to investigate the accuracy of the mode-coupling approx- 
imation by comparing it with the recently obtained exact 
solution of Prahofer and Spohn [jl0| for d = 1 . The mode- 
coupling approximation is gratifyingly close to the exact 
solution, which encourages one to believe in the utility of 
the mode-coupling approach in higher dimensions where 
no exact solution is known or likely to be found. 

In a recent paper we studied the long time properties 
of the KPZ equation within the mode-coupling approxi- 
mation, and we predicted a stretched exponential decay 
for the correlation function at long times. In this pa- 
per we have found numerically the solution of the mode- 
coupling equations in d = 1 which confirms the results of 
the previous asymptotic analysis but which also reveals 
that the correlation functions decay to zero in an oscil- 
latory manner - a fact which was not revealed by our 
previous asymptotic analysis. 

Mode-coupling equations are coupled equations for the 
correlations and response function. The correlation and 
response function are defined in Fourier space by 

COcw) = (Mk,w)fc*(k,w)>, 
G(k,w) = <flfe(k,w)/^(k,w)>, 
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where (■) indicate an average over rj. In the mode- functions are the solutions of two coupled equations, 
coupling approximation, the correlation and response 

G-\k, u) = Go" 1 (k, u) + A 2 / f / (gf, [q ■ (k - q)] [q ■ k] G(k - q, W - 0)G(q, fi) , (2) 
G(k, W ) = G (k, W ) + £ | G(k, W ) | 2 / f / ^ [q • (k - q)] 2 G(k - q, W - tt)G(q, fi) , (3) 



where Go(k,o;) = {uk 2 — iu)^ 1 is the bare response func- 
tion, and G (k, w) — 2D | G(k, w) | 2 . In the strong 
coupling limit, G(k, w) and G(k, uj) take the following 
scaling forms 

G(k,u) = k- z g{x) , 
G(k,w) = fc-( 2 x+ d+z )n(x) , 

with a; = u/k z . In d = 1, the mode-coupling equations 
simplify, due to the existence of a fluctuation dissipa- 
tion theorem which relates the correlation function to 
the response function. In t and k space, the fluctuation 
dissipation theorem can be written as: 

G(k,t) = ^-e(t)C(k,t). (4) 

(We use the same notation G and G in t space and u> 
space and indicate which one we mean by the arguments). 
We choose v = 1 and D = 1 in what follows. The mode- 
coupling equations are then reduced to one single equa- 
tion that in terms of the response function in k and t 
space G(tk z ) — G(k,t) reads: 

G(t) = 1 - A 2 f da [ dsv{s)G(o ~ s) , (5) 
Jo Jo 

where r is the scaling variable tk z , 

v(s) = — [ dxG{\ l/2 + x\ z s)G(\ 1/2 - x \ z s) , (6) 
2tt Jo 

and z = 3/2. A similar analysis to the one we have done 
can be found in |yj, where similar results were found 
(compare Fig. 1 with Fig. 1 in pl|), but here we focus 
on the long time asymptotics. To make this compari- 
son we set here and in what follows A = 1. In a pre- 
vious study we argued that an asymptotic solution 
n(r) = n^(r) for r — > oo for the scaling part of the cor- 
relation function in t and k space (r = tk z , and fi{r) is 
the Fourier transform of n(xj) is given by 

n^(r) = A (Bt) 5 , (7) 

with 7 = (d — l)/2, a = 1, and 

, 9 (0)- 2 4r(4z-4) 
P2( rf -!)/ 2 r(2z-2) 2 ' [ > 

with P = A 2 /(2'V d + 3 )/ 2 ). 
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FIG. 1. Scaling function for the response function G(r). 



In d = 1, G(r) oc h(r) due to the fluctuation dissi- 
pation theorem, so that we expect the same asymptotic 
expression for G. Our numerical analysis shows, much 
to our surprise, an oscillatory behavior for the correla- 
tion function superimposed on the stretched exponential 
decay, a feature that has never been observed before. 
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FIG. 2. Oscillation and stretched exponential behavior in 
the response function G(r). The dashed line indicates a line 
of slope —1 as predicted by Eq. (7). 
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The long time behaviour can be revealed by plotting 
\G(t)\ as a function of r 2 / 3 in a linear-log scale (see Fig. 
2), where both the presence of the oscillations and the 
overall stretched exponential decay of the envelope be- 
come apparent. We do not have any simple argument to 
explain the presence of such oscillations. However, we can 
show that they are perfectly consistent with the saddle 
point analysis performed in p2[ . The same calculation 
can in fact be repeated with a complex exponential with 
the constant B a complex number. While the calculation 
leads to the same values of 7 and a, it is not now possible 
to predict the amplitude constant A. 

We next compare the results from mode-coupling with 
the result for the scaling function in d = 1 from the exact 
solution, which recently has became available |[o|. The 
result in [fiof is given in terms of a function f(w) which 
is related to our G(r) by 



1 f°° 
f(w) = - dy 
n Jo 



cos(wy)G(y z /4) . 



(9) 



The results are shown in Fig. 3, and show a reason- 
able agreement between the mode- coupling approxima- 
tion and the exact solution p3j . From the exact solu- 
tion it is also possible to numerically calculate G(r) and 
compare it with our approximate mode-coupling solution. 
The exact solution also displays the oscillatory behavior 
which we have discovered in the mode-coupling approxi- 



mation 13 1 . 

Note that the function f(w) (g"(w)/8 in jl(J) is also 
related to the truncated correlation function in real space 



C{x,t) 



dk 
ir 



(C(fc, 0) - cos(kx)C(k, t)) = xF(t/x 3 / 2 ) 

(10) 



An earlier study of the accuracy of the mode-coupling 
approximation was undertaken by Frey et al. jy]], who 
studied the magnitude of the corrections of higher order 
diagrams to the bare vertex. They found that such di- 
agrams did produce substantial corrections. It is clear 
however that such contributions are relatively unimpor- 
tant for the correlation function we have studied. 

In summary, we have presented a numerical study 
of the mode-coupling equations for the strong coupling 
regime of the KPZ equation in the long time limit in 
d = 1 . We recovered the stretched exponential relaxation 
for the correlation function predicted previously in [JE3, 
but found a superimposed oscillation. Such oscillations 
are consistent with our previous analysis, even though we 
had not anticipated them. Finally, we compared the re- 
sults from mode-coupling theory in d = 1 with the exact 
solution from Ref. ITof . 




FIG. 3. Comparison with the results in [10]: the solid line 
is our result for f{w), the dotted line is the same function 
from Ref. [10]. 
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